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PREFACE 
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personnel of Lockheed Engineering and Management Services Company, Inc., 
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The following scientists and other personnel assisted in conducting the 
research and compiling this report. L. R. Hall and R. E. Hinkle of Lockheed 
Engineering and Management Services Company, Inc., provided assistance in the 
data-processing phases. R. B. MacDonald of the National Aeronautics and 
Space Administration, Lyndon B. Johnson Space Center, provided technical 
consultation and reviews. 

Lockheed provided, without cost, its Remote Sensing Applications Laboratory, 
located in Houston, for data processing. 
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1. SUMMARY 


The study which Is the subject of this report was conducted to explore the 
applicability of Landsat renote sensing technology to the detection of wheat- 
leaf rust In Sonora* Mexico. As a result of this study* further discussions 
with scientists of the International Center for the Improvement of Maize and 
Wheat (CIMMYT) regarding tMs problem should prove productive. 

Sonora Is Mexico's largest wheat-producing state; and Its wheat region* which 
Is located near Ciudad Obregon, Is often referred to as the "Valley of Wheat." 
Landsat data covering che Valley of Wheat during crop years 1975-76 and 
1976-77 were processed and clustered, hardcopy maps were produced* and a 
limited data analysis study was conducted In order to Identify stressed 
wheat. Crop year 1976-77 was severely affected by an epidemic of wheat-leaf 
rust* and It was believed that optimal detection of the affected wheat was 
during the 1977 jointing stage (March time frame). 

Landsat data acquired during March 1976 and March 1977 were processed, clus- 
tered* and Initially analyzed using Lockheed's Remote Sensing Applications 
Laboratory (RSAL). The Landsat multlspectral scanner (MSS) data and the RSAL 
cluster and classification results were converted Into hardcopy color maps. 
These map products were further analyzed using additional Landsat Imagery and 
cluster results (cluster means) obtained during the data- processing phases. 

Analysis of the 1977 Landsat processing results (Image maps* cluster maps, and 
cluster means) Indicated the presence of potentially stressed vegetation 
(l.e., rust-infected whtiat). Clusters containing stressed vegetation account 
for 45.6 percent of all vegetation classified as wheat using the 1977 data. 
This was significantly greater than the 13.2 percent classified using the 1976 
data. The greatest concentration of wheat classified as being stressed with 
wheat-leaf rust using the 1977 data was centered near Huatabampo* Sonora* 
which was the early center of the epidemic of wheat-leaf rust. 
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Analysis revealed a significant Increase In fallow fields In the Sonora 
agricultural region from March 1976 to March 1977. Fallow fields Increased 
from approximately 2 percent of the cropland In 1976 to approximately 
25 percent In 1977. Analysis of the data also Indicated a reduction In wheat 
area of 52 461 hectares (129 635 acres) between 1976 and 1977. 

Landsat full-frame analysis revealed that the two major reservoirs serving the 
Valley of Wheat were extremely low In 1977 when compared to water levels 
detected In 1973 Landsat Imagery. Based on this Information, It was hypothe- 
sized that the reduction In reservoir Irrigation waters and the Issues of land 
tenure disrupted the wheat and summer crop (cotton and soybean) plantings for 
the 1976-77 crop cycle, resulting In an Increase of fallow (unplanted) fields 
of 52 461 hectares (129 635 acres). A literature search of the U.S. Depart- 
ment of Agriculture (USDA) Foreign Agriculture periodicals for the years 
1975-80 provided Information to substantiate these conclusions (refs. 1-3). 

If follow-on studies are to be conducted, It Is recommended that (1) additional 
Landsat data covering the 1975-76 and 1976-77 growing cycles be collected, 

(2) Improved crop calendars be obtained, (3) some of the newly developed clas- 
sification techniques be utilized, and (4) Sonora ground-observations data 
collected by CIWYT personnel be obtained for additional analysis. 
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2. INTRODUCTION AND OBJECTIVES 


The wheat-growing region In the State of Sonora , Mexico, Mas severely affected 
by an epidemic of wheat-leaf rust during crop year 1976-77. The wheat-leaf 
rust, which Is caused by a fungus ( Puooinia rtoondita ), significantly reduced 
the state's wheat production. 

Sonora Is Mexico's largest wheat- producing state; and Its wheat region, which 
Is located near Ciudad Obregon, Is often referred to as the "Valley of Wheat." 

CIMMYT Informed the National Aeronautics and Space Administration (NASA), 
Lyndon B. Johnson Space Center (JSC), that the early center of the wheat-leaf- 
rust epidemic was at Huatabampo, Sonora, In February 1977. This town Is situ- 
ated to the southeast of Ciudad Obregon In the Valley of Wheat. The wheat- 
leaf rust then spread from Huatabampo Into the Ciudad Obregon area. By March 
when wheat was In the jointing stage, the wheat-leaf rust had affected almost 
all areas In the valley. During the March time frame, the rust was clearly 
visible and had reached a maximum in terms of external plant appearance. 

Apart from the problem of wheat- leaf rust, It was believed that Issues of land 
tenure (property ownership) In Sonora had Impacted the 1976-77 crop production 
cycle. Large commercial farms were occupied by farm laborers during 1976; and 
late In 1976, the Mexican Government declared that ownership of nearly 
100 000 hectares (247 105 acres) In Sonora was Illegal. The property rights 
against some of this land remain In dispute, and the resolution of the Issue 
Is contingent upon the decision of the Mexican Supreme Court. 

An additional and Important factor that affected crop production In Sonora was 
a drought which began In 1975. Image analysis revealed that the water levels 
In local reservoirs were extremely low In 1977 when compared to the 1973 
levels. Information from the USDA Foreign Agriculture periodicals showed that 
reduced reservoir levels Impacted crop plantings In the 1976-77 crop year. 
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In response to CIMMYT's Interest In the application of remotely sensed tandsat 
data to agricultural problems, the NASA/ JSC Space and Life Sciences Directorate 
Initiated a project to explore the use of recently developed aerospace remote 
sensing technology to the specific problem of detecting wheat-leaf rust. 

The objective of the project was to explore the applicability of machine- 
classified Landsat data to the detection of wheat-leaf rust. The specific 
objectives of the study were as follows: 

a. To cluster and classify Landsat data acquired over Sonora during March 
1977 and March 1976 

b. To produce hardcopy products from the Landsat MSS data (Landsat Image map) 
and the classified data (cluster map) 

c. To perform a limited comparative analysis of the Landsat Image maps and 
cluster maps 
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3. SONORA STUDY AREA 


The Sonora exploratory study site was limited to an area of approximately 
7770 square kilometers (3000 square miles) In the Valley of Wheat near Ciudad 
Obregon, Sonora (fig. 3-1). This represents approximately DC percent of the 
agricultural land In the Valley of Wheat. The remaining 10 percent of the 
agricultural land Is located within an adjacent Landsat footprint (37-41). A 
computer-compatible tape (CCT) set for this footprint was on order from the 
Earth Resources Observation Systems (EROS) Data Center, Sioux Falls, South 
Dakota, but was unavailable at the time data processing and analysis were 
conducted. 



Figure 3-1.- Sonora, Mexico, study area location. 
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4. SONORA AGRICULTURAL REGION 


The Sonora agricultural region Is based, largely, upon an Irrigated multiple- 
cropping system which produces a variety of fields crops (e.g., wheat, 
soybeans, and cotton), vegetables, and fruits. 

Sonora is the largest wheat -producing state in Mexico. The majority of 
Sonora's wheat and other agricultural crops are grown on the southern coastal 
plain. The northern extent of the agricultural region is near Hermosillo and 
ranges south along the Gulf of California to the State of Sinaloa. 

The coastal plain, which is generally warm and dry, is in year-round produc- 
tion of irrigated crops such as wheat and soybeans. Yields per planted unit 
rank among the highest in the world. Wheat has typically been produced by the 
more progressive farmers during the cool, dry winter months when water for 
irrigation is available (ref. 4). The source of water is runoff from the 
coastal Sierria Madre highlands during the July-October time period. The 
water is stored in large reservoirs near Ciudad Obregon. 

The study site selected for this demonstration of Landsat remote sensing tech- 
nology is situated near Ciudad Obregon. This area represents Mexico's major 
wheat region and is often referred to as the Valley of Wheat. 

The principal crops grown in Sonora are wheat, soybeans, safflower, cotton, 
sorghum, sesame, alfalfa, and flax. Summer crops such as cotton and soybeans 
are doubled-cropped with wheat. 

Cotton and soybean production in Sonora has varied greatly because of world 
cotton prices and government price supports. When cotton prices are depressed, 
the farmers have turned to soybeans. In 1977, the cotton area was supposed to 
increase in response to (1) higher world prices and (2) an expected decrease in 
soybean plantings because of reduced irrigation waters. Cotton can replace 
soybeans since its water requirements are less, thus requiring less irrigation. 
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5. WHEAT-LEAF RUST 


Wheat-leaf rust Is caused by a fungus ( Puooinia reoondita) which Is found all 
over the world (ref. 5). The fungus attacks leaf blades and sheaths, resulting 
In orange or orange-brown pustules on the green leaves. At times, the fungus 
also affects the stem, glumes, and awns. The fungus depletes the plant's water 
and nutrients needed for development of the kernel. As a consequence, the 
plant transpires water at an accelerated rate, resulting In stress. 

The red-spore stage, which is responsible for the plant damage, will result In 
fewer and lighter kernels. Yield losses frou wheat-leaf rust may be 30 per- 
cent, but losses up to 94 percent have been reported for susceptible wheat 
varieties. 

NASA/JSC was informed by CIMMYT that the use of March Landsat acquisitions 
would probably be most suitable for detecting rust-infected wheat. * 


Personal communication to R. B. MacDonald, NASA/JSC, from R. Glenn Anderson, 
director of the Wheat Program for CIMMYT, dated July 5, 1979. 
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DATA SET 


The Landsat digital data used In the exploratory study consisted of data from 
1976 and 1977. Eight CCT sets which cover four Landsat footprints were 
ordered for both 1976 and 1977. The EROS Data Center scene numbers and the 
acquisition dates for these data are as follows: 


Landsat scene 
Identification 

Acquisition 

date 

Landsat 
f ootpri nt 

8532516444500 

3/09/76 

37-41 

8242017021500 

3/17/76 

36-41 

8278116555500 

3/13/77 

37-41 

8278016501500 

3/12/77 

36-41 

8292516483500 

8/04/77 

37-41 

8292416425500 

8/03/77 

37-41 

8539716395500 

5/20/76 

37-41 

8260016570500 

9/13/76 

36-41 


The March 17, 1976, and March 12, 1977, Landsat scenes were selected for 
computer-assisted clustering and classification. The remainder of the CCT's, 
along with five full-frame color composites from 1977 and a black-and-white 
band 5 Image from 1976, were used to Identify general land uses, wheat, and 
other agricultural crops. 

The EROS Data Center scene numbers and the acquisition dates for the full- 


frame color composites are as 

follows: 


Landsat scene 

Acquisition 

Landsat 

identification 

date 

footprint 

8274516572500 

2/05/77 

37-41 

8276216510500 

2/22/77 

36-41 

8278016501500 

3/12/77 

36-41 

8278116555500 

3/13/77 

37-41 

8292516483500 

8/04/77 

37-41 
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In addition, the 1976 black-and-white Image, dated September 13, along with 
the March 17, 1976, color Landsat Image map, was used to Identify specific 
agricultural crops (e.g., wheat). The EROS Data Center scene number for the 
1976 black-and-white image is 8260016570500 (Landsat footprint 36-41). 

Multitemporal coverage of the entire study area was available only for the 
sldelap area between the footprints 36-41 and 37-41. As a result, postharvest 
data were not available for the eastern portions of the study area. The 
effects from this are discussed in section 10, Data Analysis Results. 
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7. TECHNICAL APPROACH 


The technical processing and analysis methodologies for this study are 
described In this section. A generalized flow diagram for the data processing 
and analysis functions Is shown In figure 7-1. 

7.1 REMOTE SENSING APPLICATIONS LABORATORY 

The data- processing system used to cluster and classify the Landsat MSS data 
was the Lockheed RSAL. The system Is a user-oriented digital Image-processing 
system consisting of the Digital Equipment Corporation PDP 11/70 central 
processing unit, RP06 disk pack subsystem, TU45 magnetic tape subsystem, 

LA36 Decwriter II, CR11 card reader, VT52 cathode-ray tube (CRT) terminal. 
Model 4012 graphics terminal; a Floating Point Systems, Inc., AP-1208 array 
processor; and a COMTAL Corporation Vision One Interactive image analysis 
system. 

7.2 DATA PROCESSING 

The Image analysis and computer-aided classification techniques utilized to 
produce the digital Landsat image maps and cluster maps are described in the 
following sections. 

7.2.1 INITIAL IMAGE ANALYSIS 

The Landsat full-frame Imagery of Sonora for 1973, which was available at the 
NASA/JSC Earth Observations Division, was analyzed In order to make a final 
selection of the exploratory study area. Analysis of the data revealed that 
the Valley of Wheat near Ciudad Obregon was covered In its entirety by two 
adjacent full frames (Landsat footprints 36-41 and 37-41). However, data 
availability constraints resulted in the selection of only Landsat data cover- 
ing footprint 36-41 for the computer-aided clustering and classification 
process. 
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Figure 7-1.- Generalized technical processing and analysis flow. 
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7.2.2 COMPUTER- ASSISTED CLASSIFICATION 
7. 2. 2.1 Data File Creation 


The Landsat CCT sets were loaded, sampled, and displayed on the RSAL digital 
Image display system. Data files were created from areas of approximately 
1600 columns by 1600 scan lines from the March 12, 1977, and'March 17, 1976, 
CCT's. These data were then used for detailed Image analyses and for cluster- 
ing and classification. This resulted In an Image size of approximately 25 by 
22 Inches when displayed at a scale of 1:250 000 (the scale specified by the 
contract job order). 

7. 2. 2. 2 Clustering 

For more efficient processing, the data sets were reduced by sampling the 
files for each Landsat subscene (March 1977 and March 1976). The data were 
clustered using the Iterative Self-Organizing Clustering System (ISOCLS) 
program. The RSAL ISOCLS algorithm included modifications to the original 
algorithm to reflect new requirements from a variety of users. 

An Initial unsupervised clustering of the 1976 data was conducted using an 
established maximum of 15 clusters; however, analysis of the results Indicated 
that too many classes of land use were contained in each cluster. The clus- 
tering parameters were then set to generate a maximum of 60 clusters for each 
Landsat subscene using the unsupervised mode. Clustering of the March 12, 
1977, data resulted In 43 clusters; the March 17, 1976, data clustered Into 
44 clusters. 

Following the clustering of the two sampled subscene data files, the picture 
elements (pixels) which were not sampled and clustered were then assigned to 
the appropriate cluster via the nearest-neighbor criterion. The Landsat MSS 
Image data and cluster results were then written Into the local Universal 
format, which Is compatible with the digital tape-to-fllm conversion system. 
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7. 2. 2. 3 Manual Classification Technique 


The clustering results for each crop year, 1976 and 1977, were analyzed using 
the RSAL Image analysis station. Each cluster was merged (overlaid) onto the 
Landsat MSS color-composite Image, and the land use or land uses contained 
within the clusters were determined. This procedure was used to analyze each 
cluster for both 1976 and 1977. 

Clusters containing similar land uses were noted, and lists of the 1976 and 
1977 clusters which should be aggregated and coded under a single category 
(e.g., water, woodland, grassland, and barren) were transmitted to Selscom 
Delta Inc. (Houston) for conversion Into color-coded cluster maps (see 
tables 7-1 and 7-2). 

Agricultural (cropland) land uses, especially those having a vegetative cover, 
were coded separately. The data analysis was to concentrate on detection of 
rust-infected wheat, and any cropland cluster was a candidate for containing 
the stressed wheat. Consequently, cropland clusters were coded after the 
products were generated. 
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TABLE 7-1.- CLUSTER AGGREGATIONS FOR 1976 DATA 


Cluster 

nap 

number 

Cluster name 

Original 
cluster numbers 

Noncropland clusters 

1 

Water and shadows 

8, 12, 29, 32, 39, 46 

2 

Barren land and urban 

1, 9, 16, 19, 23, 25, 27, 35 

3 

Woodl and 

4, 14, 20 

4 

Desert shrub and fallow 

10, 31 

5 

Grassland 

5, 13, 30, 34, 37, 38 

6 

Roads and border pixels 

2, 17, 21, 22, 26, 33, 36, 41, 42, 
45, 47 

Cropland clusters 

7 

Wheat 

3 

8 

Wheat 

6 

9 

Wheat 

7 

10 

Wheat and alfalfa 

11 

11 

Wheat 

15 

12 

Wheat and alfalfa 

18 

13 

Wheat 

24 

14 

Wheat borders and roads 3 

28 

15 

Wheat 

40 

16 

Wheat 

43 

17 

Wheat borders 

44 


Cropland and noncropland spectral confusion. 
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TABLE 7-2.- CLUSTER AGGREGATIONS FOR 1977 DATA 


Cluster 

map 

number 

Cluster name 

Original 
cluster numbers 

Noncropland clusters 

1 

Water and shadows 

4, 8, 12, 41 

2 

Barren land 

1, 5, 24, 29, 36, 39 

3 

Woodland 

6, 34 

4 

Desert shrub and fallow 

10, 22, 33 

5 

Grassland and fallow 

9, 18, 23, 31 

6 

Roads, border pixels, 
and streams 

20, 21, 27, 32, 35, 37, 38, 44 

Cropland clusters 

7 

Fallow, grassland, and shrub 3 

2, 13, 14, 26, 28 

8 

Wheat 

3 

9 

Fallow and residue 

5 

10 

Wheat 

7 

11 

Alfalfa and wheat 

11 

12 

Alfalfa 

16 

13 

Wheat 

19 

14 

Cropland borders 

25 

15 

Wheat and alfalfa 

30 

16 

Wheat and unknown spring 
crop (late wheat) 

40 

17 

Cropland borders 

42 

18 

Wheat and alfalfa 

43 

1 


Cropland and noncropland spectral confusion. 
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8. OUTPUT PRODUCTS 


A description of the map products and the techniques used to generate the 
final products are presented In this section. In addition, the output from 
the RSAL processing Is discussed. 

The Sonora Landsat digital data, both MSS and cluster, were converted from the 
digital tape format Into final map products by Selscom Delta Inc. A functional 
flow of the Selscom t.ape-to-fllm conversion system Is shown In figure 8-1. * 

8.1 LANDSAT IMAGE HAP 

The Landsat Image maps for 1976 and 1977 were generated using MSS bands 4 
(visible green), 5 (visible red), and 6 (near Infrared). Bands 4, 5, and 6 
were assigned blue, green, and red, respectively. In order to simulate false- 
color- Infrared Imagery. 

Band 7, which Is normally used to produce color composites, exhibited exces- 
sive data dropouts when viewed on the RSAL CRT terminal. Band 6 did not have 
data dropouts; therefore. It was selected over band 7. 

Five copies each of the Landsat Image maps for 1976 and 1977 were produced at 
the scale of 1:250 000 (with a true-north orientation of the Image). The 
requirements of the contract job order were exceeded by providing the follow- 
ing additional map annotations: bar scale, latitude-longitude grid, locator 
map, and north arrow. 

Selscom converted the data from Universal format to Integer format and applied 
geometric corrections (pixel deskewing and pixel squaring) to each Landsat sub- 
scene Image file. Because of orbital track and satellite attitude differences 
between the 1976 and 1977 full frames, the geographic positions of many ground 
features were significantly different with respect to the latitude-longitude 

7 

Tape-to-F11m Conversion System, personal communication with Brian T. Fine, 

marketing manager, Selscom Delta Inc. 
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grid. Lockheed requested that Selscom apply a rubber- sheet correction to the 
1976 data In order to provide a better fit between the two Landsat Images and 
the grid. 

Spectral Image enhancements, which include histogram evaluations and linear 
contrast stretching, were applied In order to optimize the vegetation 
signatures. 

The enhanced Integer format data were converted Into a dot matrix and filmed 
using a laser drum plotter. Color separations for each spectral band were 
produced for the 1976 and 1977 Landsat Images. 

The latitude-longitude grid was generated with a Cal coop plotter and converted 
Into a black-and-white separation (with the grid lines displayed In black). 

The map annotations (e.g., name, date, bands, and north arrow) were filmed 
using the drum plotter and spliced Into the grid overlay separation. This 
formed a single black-and-white separation which was added to the color 
separations to produce the final Landsat Image map. 

8.2 CLUSTER HAP 

The clustered data were converted from the digital format Into cluster maps 
using a procedure similar to that used to generate the Landsat Image maps. 

The primary difference In the production of the cluster maps Is that the 
clusters are contained In a single band (l.e., each cluster Is assigned to a 
digital value ranging from 0 to 256); therefore, each cluster corresponds to a 
film density with a unique color assigned to each density level. 

The cluster maps were scaled to the Landsat Image maps and given the same 
annotations with one exception; l.e., the cluster maps have color-coded chips 
along the right margin which correspond to the data clusters displayed on the 
map. 

Five copies each of the two cluster maps (1976 and 1977) were generated. 
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8.3 CLUSTERING TABULATIONS 


In addition to producing the required map products, tabulations of the clus- 
tering results (mean, standard deviations, distance between clusters, cluster 
populations, and clustering results after each ISOCLS Iteration) were 
generated using the Lockheed RSAL. 
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9. DATA ANALYSIS TECHNIQUES 


Multi temporal image analysis techniques Mere utilized to identify the land 
uses. The Landsat full-frame imagery available for analysis is presented in 
section 6, Data Set. Historical crop calendars of state-level planting and 
harvesting seasons (ref. 6) were used to identify specific crops or crop 
classes (i.e., summer crops). The study area is located in southern Sonora 
near the border of the State of Sinaloa; consequently, an interpolative crop 
calendar was developed for this transition area between the two states. The 
Sonora crop calendars are shown in figure 9-1. 

Analysis constraints included the lack of detailed historical crop calendars, 
insufficient Landsat data (i.e., no wheat planting or harvesting data), and no 
ground-observations data from previous years. 

Cluster analysis of the 1976 and 1977 cluster maps was conducted using (1) the 
comparative image analysis technique, whereby each cluster was compared to the 
corresponding areas in the 1:250 000 Landsat image map and the Landsat full- 
frame imagery, and (2) results from an evaluation of each Landsat spectral 
band cluster means, standard deviations for each band, and distances between 
clusters. 

The imagery and cluster means were analyzed with the expectation that jointed 
rust-infected wheat, if detectable, would exhibit a stressed signature. 
Stressed vegetation usually has a less vigorous signature in the near-infrared 
region than does healthy vegetation. 

A plot of spring wheat during the jointing and ripening stages is shown in 
figure 9-2. The figure illustrates that stressed wheat, which in this case is 
ripe (senescence) wheat, can be expected to have reduced reflectance and 
emittance in the near- Infrared bands (bands 6 and 7). 

Therefore, analysis for the Identification of stressed wheat was conducted 
with the premise that a reduced col or- Infrared response or low cluster mean 


9-1 
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Figure 9-1.- Historical crop calendars for Sonora 
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Figure 9-2.- Reflectance of spring wheat at different maturity 
stages for Williams County, South Dakota, In 1975. 


for bands 6 and 7 during the jointing stage Mould be an Indication of rust- 
infected wheat. Reduced color-infrared response In this context can be 
roughly equated to the lower response exhibited by the ripened spring wheat 
shown In figure 9-2. 
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10. DATA ANALYSIS RESULTS 


The Imagery and cluster analysis results are presented In this section, along 
with observations regarding the spectral data behavior, agricultural cropping 
practices, and additional factors which may have affected cropping practices 
In 1976 and 1977. 

10.1 IDENTIFICATION OF CROPS 

The crop Identification methodology (decision logic) Is presented In this sec- 
tion. The methodology for 1977 Is discussed first, followed by a discussion 
of the 1976 crop identification techniques. 

10.1.1 ANALYSIS OF 1977 DATA 

Agricultural fields (clusters) containing a vegetative cover in February and 
March with an absence of a vegetative cover (postharvest) In August were iden- 
tified as wheat. For a comparison of crop calendars and acquisitions, see 
figure 9-1. Some agricultural fields contained a ground cover in February, 
March, and August and are potentially either wheat double-cropped with a sum- 
mer crop or alfalfa (permanent vegetative cover). The uncertainty regarding 
the label for these fields is due to the unavailability of Landsat data for 
the May-June time frame when wheat is usually harvested; therefore, the clus- 
ters with ground-cover signatures on all three acquisitions could be labeled 
as either of the above. However, it is believed that these fields are most 
likely alfalfa since it is a significant crop in Sonora (20 000 hectares, 

49 421 acres, in 1975). 

One cluster contained fields which were nonvegetated in February, were vege- 
tated in March, and were again nonvegetated in August. These fields were 
labeled as a spring crop (potentially either wheat, safflower, or sorghum). 

Fields without a vegetative cover on all three acquisitions were labeled as 
fallow, although the possibility exists that a summer crop may have been sown 
following the wheat harvest In May or June and harvested prior to the August 
acquisition. 
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It should be noted that a single acquisition (March 12) was clustered and 
classified; therefore, clusters are mixed (l.e., contain more than one crop), 
as Indicated In tables 7-1 and 7-2. 

10.1.2 ANALYSIS OF 1976 DATA 

The available Landsat acquisitions for 1976 which were used In the analysis 
were acquired on March 9, March 17, and September 13. The March 17 acquisi- 
tion was used for the clustering and classification. A field with a vegeta- 
tive cover during March and an absence of a ground cover in September was 
Interpreted to be wheat, since it is the predominant crop. 

Clusters which contained fields labeled as alfalfa using the 1977 data were 
also labeled as alfalfa (wheat plus alfalfa) using the 1976 data. (Alfalfa Is 
a perennial crop and may remain In the same field as long as 7 years.) As 
Indicated above, this resulted in mixed clusters of wheat and alfalfa. Again, 
it should be noted that the alfalfa interpretation may be in error, but it Is 
believed that this is the most likely label. The September 13 acquisition 
(black-and-white full-frame Image) was somewhat useful In determining If a 
ground cover was still present in September. 

10.2 IDENTIFICATION OF WHEAT-LEAF RUST 

The image analysis criterion for identifying stressed vegetation Is the low 
col or- Infrared response. Both the 1976 and 1977 imagery exhibited signatures 
which It Is believed fit this criterion. These less-than-vigorous (stressed) 
color-infrared signatures are dark reddish brown and light red In color. 

Two 1976 wheat clusters (8 and 13) and three 1977 wheat clusters (10, 13, and 
16) have been Interpreted as containing stressed signatures. The 1977 stressed 
clusters probably contain the rust-infected wheat. The 1976 stressed clusters 
may be due to either wheat-leaf rust, a disease similar to leaf rust, or 
healthy crops with reduced Infrared reflectance and emlttance characteristics. 
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Wheat cluster means for each 1976 and 1977 landsat MSS band were analyzed, and 
the results are shown In figures 10-1 and 10-2. 3 It should be noted that the 
clusters Identified as containing the stressed-wheat signatures have the lowest 
cluster means In the near-infrared bands (6 and 7). Also, It Is Interesting to 
note that the threshold for Image Interpreted vigorous (healthy) wheat Is at 
approximately digital value 50 or less In band 7 for both crop years. The 
cluster means for noncropland are shown In figure 10-3 for comparison purposes. 

The clusters Identified as containing stressed vegetation (potentially rust- 
infected wheat) In 1977 account for 45.6 percent of the total clusters Identi- 
fied as wheat (or predominately wheat). For 1976, the stressed clusters 
account for a significantly lower proportion (only 13.2 percent) of the total 
wheat clusters. (Tables 10-1 and 10-2 contain the wheat cluster populations 
used to calculate these proportions.) 

The greatest concentration of wheat classified as being stressed with wheat- 
leaf rust using the 1977 data was centered near Huatabampo, Sonora, which was 
the early center of the epidemic of wheat-leaf rust. 

10.3 INCREASE IN FALLOW FIELDS 

A significant increase In the number of fallow fields from 1976 to 1977 was 
observed. In March 1976, approximately 2 percent of the cropland was In 
fallow; however, in March 1977, approximately 25 percent of the cropland was 
In fallow. 

10.4 REDUCTION IN PLANTED WHEAT AREA 

The area classified as wheat (including potential alfalfa area) In 1976 within 
the study area was 208 714 hectares (515 745 acres). The area sown to wheat 
in 1977 was 156 253 hectares (386 110 acres), a reduction of 52 461 hectares 
(129 635 acres). Again, the wheat area Includes potential alfalfa confusion 
area. 

3 — — 

The digital values for band 7 (near Infrared) were doubled In order to 
compensate for data compression which occurs during Landsat telemetry. 
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TABLE 10-1.- WHEAT CLUSTER POPULATIONS FOR 1976 DATA 


Cluster 

Crop 

Total pixels 
In each cluster 

Pixels In 
stress clusters 

7 

Wheat 

6 525 


8 

Wheat 

22 146 

22 146 

9 

Wheat 

12 429 


10 

Wheat and alfalfa 

114 800 


11 

Wheat 

66 054 


12 

Wheat and alfalfa 

84 140 


13 

Wheat and a sunmer crop 

40 402 

40 402 

14 

Wheat, roads, and wheat borders 

17 880 


15 

Wheat 

36 355 


16 

Wheat 

45 746 


17 

Wheat and wheat borders 



26 693 


Total pixels 

a 473 170 

b 62 548 


Represents 208 714 hectares (515 745 acres). 
Represents 27 590 hectares (68 177 acres). 
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TABLE 10-2.- WHEAT CLUSTER POPULATIONS FOR 1977 DATA 


Cluster 

Crop 

Total pixels 
In each cluster 

Pixels In 
stress clusters 

8 

Wheat 

10 714 


10 

Wheat 

43 180 

43 180 

11 

Alfalfa and wheat 

80 150 


13 

Wheat 

75 356 

75 356 

14 

Borders with crop 

18 


15 

Wheat and alfalfa 

65 185 


16 

Wheat and a spring crop 

42 994 

42 994 

17 

Borders with crop 

49 


18 

Wheat and alfalfa 

36 583 


Total pixels 

a 354 229 

b 161 530 


Represents 156 253 hectares (386 110 acres). 
Represents 71 252 hectares (176 068 acres). 
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Figure 10-1.- Plot of 1977 wheat cluster means. 
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Figure 10-2.- Plot of 1976 wheat cluster means. 
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Figure 10-3.- Plot of 1977 noncropland cluster means. 
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10.5 REDUCTION IH RESERVOIR LEVELS 

Analysis of Landsat full-frame Imagery acquired during late 1976 and 1977 Indi- 
cates reduced water levels In the reservoirs which provide Irrigation water to 
the Valley of Wheat. The water levels In late 1977. especially, are signifi- 
cantly lower than the levels detectable In the 1973 Landsat Imagery. The 
affected eservolrs are the Presa Alvara Obreyon, north of Ciudad Obregon, and 
Presa Macuzarl, east of the Valley of Wheat. 
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11. CONCLUSIONS 


Analysis of the 1977 Landsat lata Indicates that i significant proportion of 
the Identified wheat Is stressed (potentially rust-infected wheat). Addition- 
ally, a significant Increase In fallow fields from 1976 to 1977 was detected. 
Initially, It was concluded that the fallow fields represented unplanted 
cropland. It was believed that this resulted from planting disruptions due to 
land occupation by the farm workers or that it was due to a shortage of reser- 
voir Irrigation waters. The net effect from the Increased fallow fields Is a 
reduction In harvested wheat area of 52 461 hectares (129 635 acres). Subse- 
quent to the Initial Interpretation, a literature search of the USOA Foreign 
Agriculture periodicals for the years 1975-80 was conducted. Information from 
these periodicals has substantiated the conclusion regarding the effects from 
the Mexican land reformation action. It also substantiated the conclusion 
that the reduction In 1977 planted area was due, partially, to a lack of 
Irrigation waters (low reservoir levels). 

Ground observations will be required to substantiate these analyses and con- 
jectures. The possibility exists that wheat-leaf rust is not Landsat detect- 
able and that the clusters Identified as containing stressed signatures 
represent different varieties of wheat or perhaps even nonwheat crops. 

Analysis of the data was constrained by Insufficient Landsat data for fall 
1976 and summer 1977, incomplete historical crop calendars, and an absence of 
ground observations. 
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12. RECOMMENDATIONS 


it is rec ssmended that, if follow-on exploratory studies are conducted, a sub- 
set of the current study area be processed using more advanced classification 
techniques. Th*s would permit clustering and analysis of multitemporal data, 
as opposed to the single-acquisition clustering used in this study. It is 
also recommended that NASA/JSC acquire additional Landsat CCT data for the 
1975-76 and 1976- T growing cycles. The Sonora ground- observations data 
collected by ClWii' and improved crop calendars should also be acquired prior 
to performing follow-on studies. 
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